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Nonperturbative properties of gluons are studied in SU(3) lattice QCD at the quenched level. The 
first subject is a functional-form analysis of the gluon propagator D a ^ v (x) in the Landau gauge. 
We find that the gluon propagator D"^ (r) obtained in lattice QCD is well described by the four- 
dimensional (4D) Yukawa-type function e~~ mr /r with m ~ 600MeV for the Euclidean 4D distance 
r = 0.1 ~ 1.0 fm. In momentum space, the gluon propagator Zfe (p 2 ) (p — 0.5 ~ 3 GeV) is found 
to be well approximated with a new-type propagator of (p 2 + ot 2 )~ 3 / 2 , which corresponds to 4D 
Fourier image of the Yukawa-type function. Associated with the Yukawa-type gluon propagator, 
we derive analytical expressions for the zero-spatial-momentum propagator Do (f), the effective 
mass M e ff(f), and the spectral function p(») of the gluon field. The mass parameter m turns 
out to be the infrared effective mass of gluons. The obtained gluon spectral function p(ffl) is 
almost negative-definite for co > m, except for a positive 5-functional peak at 0) = m. The second 
subject is a lattice-QCD determination of the relevant gluonic momentum-component for color 
confinement. As a remarkable fact, the string tension is found to be almost unchanged even after 
cutting off the high-momentum gluon component above 1.5 GeV in the Landau gauge. In fact, 
the relevant gluonic scale for color confinement is concluded to be below 1.5 GeV. 
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1. Introduction 

Quantum chromodynamics (QCD) and the gluon field A^(x) = A^(x)r fl G su(3) were first 
proposed by Nambu in 1966 [jl|] just after the introduction of color degrees of freedom. Although 
QCD has been established as the fundamental gauge theory of the strong interaction with many 
successes, there are still unsolved problems on nonperturbative QCD in the low-energy region. 

The analysis of gluon properties is an important key point to clarify the nonperturbative aspects 
of QCD [gp. In particular, the gluon propagator, i.e., the two-point Green function is one of the 
most basic quantities in QCD, and has been investigated with much interests in various gauges, 
such as the Landau gauge [|, § |, |, 0, § g, |Kj, [TT], 0, [B|, |§, the Coulomb gauge ||B|, ||], and 



the maximally Abelian (MA) gauge [|l7|, |18| ], in the context of various aspects of QCD. Dynamical 
gluon-mass generation [Q,[[9|] is also an important subject related to the infrared gluon propagation. 
While gluons are perturbatively massless, they are conjectured to acquire a large effective mass as 
the self-energy through their self-interaction in a nonperturbative manner. For example, glueballs, 
color-singlet bound states of gluons, are considered to be fairly massive, e.g., about 1.5GeV for the 



lowest ++ and about 2GeV for the lowest 2 ++ , as indicated in lattice QCD calculations [ ]20| , [21]]. 

In this paper, we study the functional form of the gluon propagator in the Landau gauge in 
SU(3) lattice QCD Monte Carlo calculations, especially for the infrared and intermediate region of 



r = 0.1 ~ 1 .Ofm, which is relevant for the quark-hadron physics [gOj, |22p, and aim at a nonperturba- 
tive description of gluon properties, based on the obtained function form of the gluon propagator. 
As another subject, using lattice QCD, we also study the relevant gluonic momentum-component 



for color confinement at the quantitative level, by introducing a cut in the momentum space [23]. 



2. Formalism for gluon propagator in Landau gauge 

In this section, we briefly review the formalism of Landau gauge fixing and the gluon prop- 
agator in Euclidean space-time. The Landau gauge is one of the most popular gauges in QCD, 
and keeps Lorentz covariance and global S\J(N C ) symmetry. Owing to these symmetries and the 
transverse property, the color and Lorentz structure of the gluon propagator is uniquely determined. 

In Euclidean QCD, the Landau gauge has a global definition to minimize the global quantity, 

R = J d 4 x Tr{A^ (x)A jU (x)} = - J d 4 xA a ^ (x)A£ (x) , (2.1) 

by gauge transformation. The local condition d^A^ (x) = is derived from the minimization of R. 
The global quantity R can be regarded as "total amount of the gauge-field fluctuation" in Euclidean 
space-time. In the global definition, the Landau gauge has a clear physical interpretation that it 
maximally suppresses artificial gauge-field fluctuations relating to gauge degrees of freedom [jH|]. 
We then expect that only the minimal quantum fluctuation of the gluon field survives in the Landau 
gauge, and the physical essence of gluon properties can be investigated without suffering from 
artificial fluctuations of gauge degrees of freedom. 

In lattice QCD, the gauge field is described by the link- variable Up(x) = e' agAti ^ with the 
lattice spacing a and QCD gauge coupling g. The Landau gauge is defined by the maximization of 

/?iatt = ££ReTrt/ M (x), (2.2) 

x n 
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by the gauge transformation. The maximization of R\ att corresponds to the minimization of R in the 
continuum theory, and maximally suppresses the gauge-field fluctuation. 

The bare gluon field A b * Te (x) G su(N c ) is defined from the link-variable as 

< e (*) = ^ [U^x) - Ufa)] - 2-^-Tr [U,(x) - Ufa)] , (2.3) 

where the second term is added to make A b ™ e traceless. In the Landau gauge, the minimization of 
gluon-field fluctuations justifies the expansion by small lattice spacing a. The renormalized gluon 
field A^ (x) is obtained by multiplying a real renormalization factor Z 3 ^ 2 as A^ (x) = Z 3 1 ' 2 A^ K (x) . 
In the Euclidean metric, the gluon propagator D a ^ v (x) is defined by the two-point function as 

Dfvfry) = {A a il {x)A b v (y))=Df v {x-y). (2.4) 
In the coordinate space, we investigate the scalar combination of the gluon propagator 

m - 3(jy^W = 3^1) WW®)' ^ 
as a function of the 4D Euclidean distance r = \x\ = {x^Xp) 1 ! 2 . 



3. Lattice QCD result for the functional form of gluon propagator in Landau gauge 

We study the functional form of the gluon propagator D(r) = DffJr) /24 in the Landau gauge 
in SU(3) lattice QCD, in the infrared and intermediate region of r = {x a x a ) l l 2 = 0.1 ~ l.Ofm, 
which is the relevant scale of quark-hadron physics. Here, we mainly deal with the coordinate- 
space gluon propagator D(r), which is directly obtained from lattice calculations and then more 
primary than the momentum-space propagator D(p 2 ) = J d 4 x e' px D{x). The SU(3) lattice QCD 
Monte Carlo calculations are performed at the quenched level using the standard plaquette action 
with /3 = 2N c /g 2 =5.1, 5.8, and 6.0, on the lattice size of 16 3 x 32, 20 3 x 32, and 32 4 , respectively. 
The lattice spacing a is found to be a = 0.186,0.152, and 0.104fm, at j8 = 5.7, 5.8, and 6.0, 
respectively, when the scale is determined so as to reproduce the string tension as y/a = 427MeV 



from the static QQ potential Q24fJ . Here, we choose the renormalization scale at \x = 4GeV for 
j3 = 6.0 [|| ^J, and make corresponding rescaling for j8=5.7 and 5.8. 

Figure 1(a) and (b) show the coordinate-space gluon propagator D(r) and the momentum- 
space gluon propagator D(p 2 ), respectively. Our lattice QCD result of D(p 2 ) is consistent with that 



obtained in the previous lattice studies []8|, |9|], although recent huge- volume lattice studies QllJ , |12[ ] 
indicate a suppression of the gluon propagator in the Deep-IR region (p < 0.5GeV), compared with 
the smaller lattice result. 



First, we consider the coordinate-space propagator of the free massive- vector field [17], 



, , f d 4 p ,„ _ 1 1 m 

D mass (r) = / j^e ~'P - x ^— - = —- Kl (mr), (3.1) 

with the modified Bessel function K\(mr). For large r, one finds K\{mr) ~ a/ 2^?^ ^' anc ^ tne 
massive propagator behaves as Z) mass (r) ~ r~^l 2 e~ mr . In Fig. 1(a), we add the fit result of the lattice 
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0.2 0.4 0.6 0.8 1 1 2 3 4 

r [fm] p [GeV] 



Figure 1: (a) Lattice QCD data of the Landau-gauge gluon propagator D(r) = D^ fl (x)/24 in coordinate 
space at j3=5.7, 5.8, and 6.0, and the Yukawa-type function D-ftikawaM = Ame~ mr / V (solid line) with m = 
0.624GeV and A = 0. 162. The dash-dotted line denotes a typical example of the massive-vector propagator 
D mass (r). (b) The Landau-gauge gluon propagator D(p 2 ) — S^(/? 2 )/24 in momentum space. The symbols 
denote the lattice QCD data at /3 = 6.0, where the momentum is defined as = | sin(-^). The solid line 
denotes the Yukawa-type propagator in the momentum space, i.e., £>Yukawa(/? 2 ) = An 2 Am(p 2 + m 2 )~ 3 / 2 . 



data with D mass (r) in the fit-range of r = 0.6 ~ l.Ofm. In the IR region, this fit seems well, and the 
effective mass m is estimated to be about 500MeV. However, as shown in Fig. 1(a), the lattice gluon 
propagator D(r) cannot be described with D mass (r) in the whole region of r = 0.1 ~ l.Ofm. 

By the functional-form analysis of the gluon propagator, we find that the Landau-gauge gluon 
propagator D(r) in the coordinate space is well described by the 4D Yukawa-type function [14] 

e -mr 

OYukawaO") =Am (3.2) 

r 

with m = 0.624(8)GeV and A = 0.162(2) in the range of r = 0.1 ~ l.Ofm, as shown in Fig. 1(a). 

In Fig. 1(b), we add by the solid line the Fourier transformation of the Yukawa-type function 
-DYukawafr), i.e., Dy u kawa(^ 2 ) = 4-7t 2 Am(p 2 + ra 2 )~ 3 / 2 , with the same parameters m = 0.624GeV 
and A = 0.162 as those used for the coordinate-space gluon propagator. The lattice QCD data of 
D(p 2 ) are found to be approximated with Yu kawa(/' 2 ) m the range of p < 3GeV, while they seem 
to be consistent with the tree-level massless propagator D tree (p 2 ) = l/p 2 for large p 2 . 

We summarize the functional form of the gluon propagator obtained in SU(3) lattice QCD: 

1. The gluon propagator D(r) in the Landau gauge is well described by the four-dimensional 
(4D) Yukawa-type function as [ 14 ] 

1 e~ mr 
D(r) = ^D a w(r)=Am—, (3.3) 

with m ~ 600MeV and A ~ 0.16, in the whole region of r = (x a x a ) 1//2 = 0.1 ~ l.Ofm. 

2. The gluon propagator D{p 2 ) in the momentum space is also well described by the corre- 
sponding new-type propagator (4D Fourier transformed Yukawa-type function) as [ 14 ] 

D( p 2) = Lo aa (p 2 ) = A%lAm (3 4) 

with m ~ 600MeV and A ~ 0.16, in the momentum region of 0.5GeV < p < 3GeV. 
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4. Analytical applications of Yukawa-type gluon propagator 



In this section, as applications of the Yukawa-type gluon propagator, we derive analytical 
expressions for the zero-spatial-momentum propagator Do(t), the effective mass M e s(t), and the 



spectral function p(ffl) of the gluon field [14Q. All the derivations can be analytically performed, 
starting from the Yukawa-type gluon propagator Dyukawa(^)- 

4.1 Zero-spatial-momentum propagator of gluons 

First, we consider the zero-spatial-momentum propagator D${t), associated with the Yukawa- 
type propagator DyukawaM, where r is the 4D Euclidean distance, r = Vx 2 + t 2 . We define the 
zero-spatial-momentum propagator Dq(/) of gluons as 



D o(0 = (0,0)) =£D(r), 



(4.1) 



where the spatial momentum is projected to be zero. For the simple argument, we here deal with the 
continuum formalism with infinite space-time. Starting from the Yukawa-type gluon propagator, 



Am Am 

^Yukawa (?) = e = < 

r y/W+t 1 

we derive the zero-spatial-momentum propagator as [ |I4| ] 



(4.2) 



D o (0 = J J 3 xD Yukawa (r) =47lAm^ 



dx 



-m\Zx 2j rt 2 



■ 4%AtK\ (mt) 



(4.3) 



In Fig.2(a), we show the theoretical curve of Do(t) in Eq.(i_3) with m=0.624GeV and A=0.162, 
together with the lattice QCD result of A)(0 in the Landau gauge. For the actual comparison with 



the lattice data, we take account of the temporal periodicity [14]. The lattice QCD data are found 
to be well described by the theoretical curve, associated with the Yukawa-type gluon propagator. 




0.8 



0.2 



M eff (t), P = 6.0 
theoretical curve 




0.4 0.6 
t[fm] 



0.8 



Figure 2: (a) The zero-spatial-momentum propagator Do (f ) of gluons in the Landau gauge, (b) The effective 
mass M e ff (f ) of gluons in the Landau gauge. The symbols are the lattice QCD data on 32 4 at j3 = 6.0, and the 
solid line is the theoretical curve derived from the Yukawa-type propagator with w=0.624GeV and A=0.162. 
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4.2 Effective mass of gluons 

Second, we investigate the effective mass M e ff(f) of gluons. This method is often used for 



hadrons as a standard mass measurement in lattice QCD [20]. For the simple notation, we use the 
lattice unit of a = 1 in this subsection. In the case of large temporal lattice size, the effective mass 
of gluons is defined as 

Meff(f) = ln{Do(/)/A)(/ + !)}• (4.4) 

In Fig.2(b), we show the lattice result of M e ff(/), where we take account of the temporal peri- 
odicity. The effective gluon mass exhibits a significant scale-dependence, and it takes a small value 
at short distances. Quantitatively, the effective gluon mass is estimated to be about 400 ~ 600MeV 
in the infrared region of about lfm p4)|. This value seems consistent with the gluon mass suggested 
by Cornwall [^], from a systematic analysis of nonperturbative QCD phenomena. 

Now, we consider the consequence of the Yukawa-type propagator DYukawa(r). For simplicity, 
we here treat the three-dimensional space as a continuous infinite- volume space, while the temporal 
variable t is discrete. When the temporal periodicity can be neglected, we obtain an analytical 



expression of the effective mass [14], 



D (t) tKAmt) 

M ef{ (t) = In — P^— = In n , -. (4.5) 

ettW A)(f + 1) (f + l)Zi(m(f + l)) 



In Fig. 2(b), we add by the solid line the theoretical curve of M^(t) in Eq.( |4.5| ) with m=0.624GeV. 
The lattice QCD data of M e s(t) are found to be well described by the theoretical curve derived 
from the Yukawa-type gluon propagator. 

From the asymptotic form K[(z) °< z7 x l 2 e~ z , the effective mass of gluons is approximated as 

M e ff(?) ~m— -ln(l + -) ~m— — (4.6) 



for large t [14]. This functional form indicates that M e ff(?) is an increasing function and approaches 
m from below, as t increases. Then, the mass parameter m ~ 600MeV in the Yukawa-type gluon 
propagator has a definite physical meaning of the effective gluon mass in the infrared region. 

Note that the simple analytical expression reproduces the anomalous increasing behavior of 
the effective mass M e ff(?) of gluons. Thus, this framework with the Yukawa-type gluon propagator 
gives an analytical and quantitative method, and is found to well reproduce the lattice QCD result. 

4.3 Spectral function of gluons in the Landau gauge 

As a general argument, an increasing behavior of the effective mass M e s(t) means that the 
spectral function is not positive-definite 0,0,0]. More precisely, the increasing property of M e g(t) 
can be realized, only when there is some suitable coexistence of positive- and negative-value re- 
gions in the spectral function p(co) [|i4[]. However, the functional form of the spectral function of 
the gluon field is not yet known. 

From the analytical expression of the zero-spatial-momentum propagator Do(t) = 4nAtK\ (mt), 
we can derive the spectral function p(ffl) of the gluon field, associated with the Yukawa-type gluon 



propagator [14]. For simplicity, we take continuum formalism with infinite space-time. 
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The relation between the spectral function p(co) and the zero-spatial-momentum propagator 
Do(?) is given by the Laplace transformation, 



D (t) 



dco p(ffl) e 



(4.7) 



When the spectral function is given by a 8 -function such asp(ft))~5(ft) — (Bo), which corresponds 
to a single mass spectrum, one finds a familiar relation of Do(t) ~ e _£B °'. For the physical state, the 
spectral function p (to) gives a probability factor, and is non-negative definite in the whole region 
of ft). This property is related to the unitarity of the S-matrix. 

From an integral representation of the modified Bessel function, we derive the following for- 
mulae on the inverse Laplace transformation, 



1 fc+i°° r,) 



1 



1 



1 



(4.8) 



2sL. * e " m,) = -w^w- e{<0 ~ 1 - e) + w^w S[m - 1 " e) - <49) 

where an infinitesimal positive e is introduced for a regularization [ITl]. Then, starting from the 



Yukawa-type propagator, we derive the spectral function p(ffl) of the gluon field as [14] 

AnAm _„ , AnA/\flm 



p{co) 



-Q(a — m — e) + - 



8(a> — m — e), 



(4.10) 



(ft) 2 — m 2 ) 3 / 2 (ft) — m) 1 / 2 

with an infinitesimal positive e. Here, m ~ 600MeV is the mass parameter in the Yukawa-type 
function for the Landau-gauge gluon propagator. The first term expresses a negative continuum 
spectrum, and the second term a 5-functional peak with the residue including a positive infinite 
factor as £ -1 / 2 at ft) = m + £. 







i i i 
p(co) 


1 






^ 

Jl 

Q. 










-2 













1 2 3 4 5 

CO [m] 

Figure 3: The spectral function p (©) of the gluon field, associated with the Yukawa-type propagator. The 
unit is normalized by the mass parameter m ~ 600MeV. p(co) shows anomalous behaviors: it has a positive 
5-functional peak with the residue of +°° at co = m, and takes negative values for all the region of CO > m. 



We show in Fig.3 the spectral function p(ft)) of the gluon field. As a remarkable fact, the 
obtained gluon spectral function p(ft)) is negative-definite for all the region of ft) > m, except 
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for the positive 5-functional peak at ft) = m. The negative property of the spectral function in 
coexistence with the positive peak leads to the anomalous increasing behavior of the effective mass 
M e ff(f) of gluons [14]. Actually, the resulting effective mass M e ff(?) well describes the lattice result, 
as shown in Fig.2(b). 

We note that the gluon spectral function p(to) is divergent at ft) = m + e, and the divergence 
structure consists of two ingredients: a 5-functional peak with a positive infinite residue and a neg- 
ative wider power-damping peak. On finite-volume lattices, these singularities are to be smeared, 
and p(ft)) is expected to take a finite value everywhere on ft). On the lattice, we conjecture that 
the spectral function p(ft)) includes a narrow positive peak stemming from the 5-function in the 
vicinity of ft) = m (+e) and a wider negative peak near ft) ~ m in the region of ft) > m [14]. 

In this way, the Yukawa-type gluon propagator indicates an extremely anomalous spectral 
function of the gluon field in the Landau gauge. The obtained gluon spectral function p(ft)) is 
negative almost everywhere, and includes a complicated divergence structure near the "anomalous 
threshold", co = m (+e). Thus, this framework with the Yukawa-type gluon propagator gives an 
analytical and concrete expression for the gluon spectral function p(ft)) at the quantitative level. 



5. A hypothesis of an effective dimensional reduction in stochastic gluonic vacuum 
by the Parisi-Sourlas mechanism 



We discuss the Yukawa-type gluon propagation and a possible dimensional reduction due to 



the stochastic behavior of the gluon field in the infrared region [14]. As shown before, the Landau- 
gauge gluon propagator is well described by the Yukawa function in /owr-dimensional Euclidean 
space-time. However, the Yukawa function e~ mr /r is a natural form in ?/jree-dimensional Euclidean 
space-time, since it is obtained by the three-dimensional Fourier transformation of the ordinary 
massive propagator (p 2 +m 2 ) -1 . In fact, the Yukawa-type propagator has a "three-dimensional" 
property. In this sense, as an interesting possibility, we propose to interpret this Yukawa-type 
behavior of the gluon propagation as an "effective reduction of the space-time dimension". 

Such a "dimensional reduction" sometimes occurs in stochastic systems, as Parisi and Sourlas 



pointed out for the spin system in a random magnetic field [|25[]. In fact, on the infrared dominant 
diagrams, the D-dimensional system coupled to the Gaussian-random external field is equivalent 
to the (D — 2) -dimensional system without the external field. 




color-magnetic fields 



propagating gluon 



Figure 4: A schematic figure for a propagating gluon in the QCD vacuum. The QCD vacuum is filled 
with color-magnetic fields which are stochastic at an infrared scale, and the gluon propagates in the random 
color-magnetic fields. 
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We note that the gluon propagation in the QCD vacuum resembles the situation of the system 
coupled to the stochastic external field. Actually, as is indicated by a large positive value of the 
gluon condensate {G a ^ v G a ^ v ) > in the Minkowski space, the QCD vacuum is filled with a strong 
color-magnetic field [|^, 26], which can contribute spontaneous chiral-symmetry breaking [|7|], 



and the color-magnetic field is considered to be highly random at the infrared scale Q26| , [280 ■ Since 
gluons interact each other, the propagating gluon is violently scattered by the other gluon fields 
which are randomly condensed in the QCD vacuum at the infrared scale, as shown in Fig.4. 

Actually at the infrared scale, the gluon field shows a strong randomness due to the strong 
interaction, and this infrared strong randomness is considered to be responsible for color confine- 
ment, as is indicated in strong-coupling lattice QCD [f2C|]. Even after the removal of the fake gauge 
degrees of freedom by gauge fixing, the gluon field exhibits a strong randomness [^] accompany- 
ing a quite large fluctuation at the infrared scale. 

As a generalization of the Parisi-Sourlas mechanism, we conjecture that the infrared structure 
of a theory in the presence of the quasi-random external field in higher-dimensional space-time has 
a similarity to the theory without the external field in lower-dimensional space-time [p~4}]. From this 
point of view, the Yukawa-type behavior of gluon propagation may indicate an "effective reduction 
of space-time dimension" by one, due to the stochastic interaction between the propagating gluon 
and the other gluon fields in the QCD vacuum, of which net physical fluctuation is highly random 
at the infrared scale. 

6. Relevant gluonic scale for color confinement 

In this section, as another subject, we study lattice-QCD analysis for the relevant gluonic 
momentum-component for color confinement Here, we formulate a new general lattice frame- 
work to extract the relevant gluonic energy scale of each QCD phenomenon by introducing a cut 



for link- variables in momentum space [23]. Our method consists of the following five steps. 



Step 1. Generation of coordinate-space link-variable in the Landau gauge 

As usual, we generate a gauge configuration on a L 4 lattice with the lattice spacing a by the 
lattice-QCD Monte Carlo simulation under space-time periodic boundary conditions, and obtain a 
finite number of coordinate-space link-variables. Here, we consider the link-variables fixed in the 
Landau gauge, which gives a transparent connection between the link- variable and the gauge field, 
owing to the suppression of gluon-field fluctuations. 

Step 2. Four-dimensional discrete Fourier transformation 

By the discrete Fourier transformation, we define the momentum-space link-variable, 

Up{p) = —— yU^{x)exTp{q^ v p v x v ), (6.1) 

"site x 

where N s i te is the total number of lattice sites. The momentum-space lattice spacing is given by 

2n 

ap = La' 
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Step 3. Imposing a cut in the momentum space 

We impose a cut on U^(p) in a certain region of the momentum space, as schematically 
shown in Fig.5. Outside the cut, we replace U^(p) by the free-field link- variable, Ojj es (p) = 
at - L.v 1 exp^XvPv^v) = 8po, corresponding to A ll (x) = or U^{x) = 1. Then, the momentum- 
space link-variable U^(p) with the cut is defined as 

U A (v) = l °^ {P) ( insidecut ) (63) 
\ujr(p) = 8 p0 (outside cut). K ' 




Figure 5: A schematic figure of the UV cut Auv and the IR cut Air on momentum-space lattice, with the 
lattice spacing a p = 2%/{Ld). The momentum-space link-variable Ufi(p) is replaced by the free variable 
t7^ ree (p) in the shaded cut regions. 



Step 4. Inverse Fourier transformation 

To return to coordinate space, we carry out the inverse Fourier transformation as 

fjto =£^(p)exp(-i£ v /7 v x v ). (6.4) 

p 

Since this U' (x) is not an SU(3) matrix, we project it onto an SU(3) element Ufr(x) by maximizing 
ReTr[L^(x)^L^,(jc)]. Such a projection is often used in lattice QCD algorithms. By this projection, 
we obtain the coordinate-space link- variable U^(x) with the cut, which is an SU(3) matrix and has 
the maximal overlap to U'^(x). 

Step 5. Calculation of physical quantities 

Using the cut link- variable Ufr(x), instead of U^{x), we calculate physical quantities as the 
expectation value in the same way as original lattice QCD. 



With this method in lattice-QCD framework, we quantitatively determine the relevant energy 
scale of color confinement, through the analyses of the QQ potential. The lattice QCD Monte Carlo 
simulations are performed on 16 4 lattice at j6=5.7, 5.8, and 6.0 at the quenched level [23]. 
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0.5 



A| R /a p = 1 i 
AWaJ= .1 i 




Figure 6: (a) The QQ potential V (R) with the IR cut Air plotted against the interquark distance R. (b) The 
QQ potential with the UV cut Auv- The lattice QCD calculation is performed on 16 4 lattice with j3 = 6.0, 
i.e., a ~ 0.10 fm and a p = 2n/ (La) ~ 0.77 GeV. The broken line is the original QQ potential in lattice QCD. 



Figure 6 (a) and (b) show the QQ potential V (R) with the IR cutoff Air and the UV cutoff Atjv, 
respectively. We get the following lattice-QCD results on the role of gluon momentum components. 

• By the IR cutoff Air, as shown in Fig. 6(a), the Coulomb potential seems to be unchanged, 



but the confinement potential is largely reduced 023Q . 



By the UV cutoff Auv, as shown in Fig.6(b), the Coulomb potential is largely reduced, but 



the confinement potential is almost unchanged [23] 



p=6.0 i— e— i 
P=5.8 i— i 

* 

l%-< 
Hi* 







0.5 1 1.5 2 2.5 3 3.5 4 
A uv [GeV] 



Figure 7: The Auv -dependence of the string tension a obtained from the asymptotic slope of the QQ 
potential V(R) with the UV cutoff Auv- The lattice QCD calculations are performed on 16 4 lattice with j3 
=5.7, 5.8 and 6.0. The vertical error-bar is the statistical error, and the horizontal error-bar the range from 
the discrete momentum. The broken line denotes the original value of the string tension a ~ 0.89 GeV/fm. 



Fig.7 shows the Auv -dependence of the string tension a obtained from the asymptotic slope of 
the QQ potential V(R) with the UV cutoff Auv- As a remarkable fact, the string tension is almost 
unchanged even after cutting off the high-momentum gluon component above 1.5 GeV. In fact, the 
relevant gluonic scale of color confinement is concluded to be below 1.5 GeV [[23|]. 
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